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Abstract. - Effective giant spins of magnetic nanoparticles are considered classically in the con- 
ventional theory of superparamagnetism based on the Landau-Lifshitz-Langevin equation. How- 
ever, microscopic calculations for a large spin with uniaxial anisotropy, coupled to the lattice via 
the simplest generic mechanism, show that the results of the conventional theory are not repro- 
duced in the limit S — > oo. In particular, the prefactor To in the Arrhenius escape rate over 
the barrier F = To exp [—At// (feT)] has an anomalously large sensitivity to symmetry-breaking 
interactions such as transverse field. 
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O Ferromagnetic particles of a sufficiently small size (e.g., 
, ^ m agnetic nanoparticles) are in a single-domain magnetic 
state, atomic spins being kept collinear by a strong ex- 
change interaction. The resulting giant spin of a mag- 
^rietic particle shows a bistability in the case of uniaxial 
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anisotropy that creates two energy minima and a barrier 
between them [1,2]. At thermal equilibrium, there is a 
distribution over directions of particles' spins similar to 
^Qthat of paramagnets. Since total spins of magnetic par- 
tides are very large, this kind of paramagnetism is called 
OO "superparamagnetism" . 

^P. Neel suggested a model of relaxation of ensembles of 
^ magnetic particles in which spins are hopping between the 
>\> tjwo energy minima [3] . Modern approach to superpara- 
; i magnetic dynamics is based on the Landau-Lifshitz equa- 
te fjion [4] for classical spins of fixed length augmented by the 
stochastic Langevin field simulating the environment [5] 
that is equivalent to the Fokker-Planck equation (FPE) 
for classical spins. Solution of the FPE yields the Arrhe- 
nius thermal activation rate T = To exp [-AU/ (fceT)] for 
T <C AU/k-B, AC/ being the energy barrier [5,6]. 

The amount of theoretical papers published on the sub- 
ject up to now is innumerable. The reader can refer to 
the book [7] on the Langevin approach to magnetic and 
dipolar systems and to Ref. [8] for a review of spin thermal 
activation problems. Numerically one can solve the FPE 
using matrix continued fractions [9] or other methods. Al- 
ternatively, one can start with the underlying stochastic 
model and solve it with matrix continued fractions [10] or 



directly as a stochastic differential equation [11,12], also 
for the model with a variable spin length near the Curie 
temperature [13]. For a model of many classical atomic 
spins forming a nanoparticlc, direct solution of a system 
of Landau-Lifshitz-Langevin (LLL) equations is the only 
working numerical method [14]. 

With increasing the size of a magnetic particle, non- 
collinearities due to the surface anisotropy proliferate and 
the single-domain state is gradually destroyed. Small 
noncollinearitics result in an additional cubic anisotropy 
for particle's global magnetization [15-17] that is second- 
order in the surfacc-anisotropy constant and scales as the 
particle's volume. Another interesting effect is relaxation 
via internal spin waves [18, 19] that should add up with 
the relaxation due to the environment. 

Experimentally, only recent successes in fabrication of 
nanoparticles with well-controlled parameters allowed to 
obtain the famous Stoner-Wohlfarth astroid and check 
the Neel-Brown theory of thermal activation [20-22] . On 
smaller nanoparticles, indications of spin tunneling [23-26] 
have been seen [27]. Later, however, the interest in spin 
tunneling has shifted to molecular magnets such as Mni2 
and Fes, where the molecular spin is only S = 10 and 
the phenomenon could be observed with a much greater 
certainty and resolution [28-32]. 

The stochastic model of magnetic particles using the 
Landau-Lifshitz equation for a large spin with the for- 
mal Langevin magnetic field has been perpetuated in the 
literature because of its simplicity. However, this model 
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contradicts the time-reversal symmetry. Deformations of 
the lattice due to thermal fluctuations cannot produce 
any fluctuating effective magnetic field (i.e., terms in the 
Hamiltonian linear in spin components). It rather pro- 
duces a fluctuating anisotropy, i.e., stochastic terms even 
in components of the spin S. The corresponding analy- 
sis has been done in Rcf. [33] where it was shown how 
the symmetry and strength of the relaxation term in 
the Landau-Lifshitz equation follows from those of the 
stochastic terms. However, this model for classical spins 
was never used because it can include too many difficult- 
to-define coupling and damping constants. 

On the other hand, all microscopic quantum-mechanical 
models of spin-lattice relaxation employ spin-lattice cou- 
plings that do not violate basic symmetries. Until re- 
cently, however, these calculations suffered from too many 
unknown coupling constants that allowed only order-of- 
magnitudc estimations. Discovery of the universal mech- 
anism of spin relaxation via distortionless rotation of the 
crystal field by transverse phonons [34,35] changed the sit- 
uation. Within this mechanism, spin-lattice coupling can 
be expressed through the parameters of the crystal field 
that can be easily measured. Implementing this mech- 
anism in the stochastic formalism for classical spins of 
Ref. [33] would allow to rewrite the theory of superpara- 
magnetism in a more satisfactory form. 

The less ambitious aim of this Letter is, however, just 
to demonstrate that relaxation of large spins of magnetic 
particles cannot be described by the conventional clas- 
sical approach. The point is that the most important 
quantum-mechanical relaxation processes such as emis- 
sion/absorption of phonons are sensitive to the energy 
levels of the spin. Since in the existing LLL formalism 
the information of the energy levels is lost, there is no 
connection to the underlying quantum mechanics and the 
ensuing results are questionable. 

One can argue that giant spins of magnetic particles, 
S ;g> 1 , are classical to a high degree of precision. This is 
not true, however, since even the relaxation in the bulk is 
governed by quantum mechanics. Of course, equilibrium 
properties of superparamagnets are classical since one has 
the Langevin function instead of the Brillouin function 
for the field-dependent magnetization. The relaxation re- 
mains non-classical, however, whatever large is the par- 
ticle. Indeed, quantum effects in magnetic particles have 
been recently observed and discussed in Ref. [36]. 

To understand the importance of quantum effects in 
magnetic particles, one has to realize the difference be- 
tween the classical- spin limit and the large-spin limit. The 
classical-spin limit is a theoretical trick to simplify calcula- 
tions by eliminating quantum effects. The large-spin limit, 
to the contrary, is the real situation. 

For instance, for a system of N atomic spins s with 
the easy-axis Hamiltonian H = —ds 2 , held together by a 
strong exchange, the total spin is S = Ns 3S> I. Within 
the classical-spin limit, the effective classical Hamiltonian 
of the system would be H = -DS 2 with D = d/N 2 



to preserve the energy barrier AU = ds 2 = DS 2 . With 
the energy levels of the spin S given by e m = —Dm 2 , 
the transition frequency ws-i,s = £s-i — £s between the 
ground and first excited states of the effective spin be- 
comes hu>s,s-i = (25 — 1) D = 2sd/N, disappearing in 
the limit N — > oo. Accordingly, the direct phonon pro- 
cesses die out for large N, so that the only relaxation 
processes due to phonons become the two-phonon Raman 
processes that become insensitive to the energy levels for 
small transition frequencies. 

In the realistic large-spin limit, the transition frequency 
u>s,s—i is preserved, because this is the frequency of the 
small-amplitude spin precession in the anisotropy field. 
Thus the effective anisotropy constant D scales accord- 
ing to huj s ,s-i = (25 - l)D = 2SD = 2NsD = 2sd = 
const, hence D = d/N. In this case the energy barrier is 
AU = DS 2 = Nds 2 cc N, the size of the particle, as it 
should be. One can see that direct spin-phonon processes 
survive in the large-spin limit. This makes the situation 
completely different from the classical-spin limit, regard- 
ing the relaxation. 

On the other hand, the transition frequencies between 
the levels near the top of the barrier, m ~ 1, are of order 
f^m.m-i = (2m — I) D ~ D oc 1/iV and they vanish in 
the large-spin limit. This means that direct phonon pro- 
cesses between the adjacent evergy levels, having the rate 

^m.m-i w m,m-i f° r ^m,m-i < k B T, die out near the 
top of the barrier that becomes a bottleneck for the ther- 
mal activation process. In this region, diffusion of spin 
populations over the stairway of adjacent levels is effec- 
tuated by much weaker Raman processes that leads to 
small escape prcfactors Tq with essential temperature de- 
pendence [37]. 

Transverse magnetic field H± or transverse anisotropy 
create saddles in the potential landscape of the effective 
spin that strongly change dynamics of thermal activation. 
A "phase diagram" of different regimes, such as uniaxial, 
high-, intermediate-, and low-damping regimes, created 
by the transverse field, has been obtained in Ref. [38]. 
Especially in the low-damping (LD) case, transverse field 
results in a strong increase of the escape rate T. As can 
be seen from the comparison of the LD and HD cases in 
Fig. 3 of Ref. [38], the main effect is the increase of the 
prefactor Tq, while lowering the barrier AU (equal in the 
LD and HD cases) plays a secondary role. 

For a quantum large spin, the effect of transverse field 
H± should be even greater, since for H± = the pref- 
actor Tq is anomalously small. For H± ^ 0, the states 
|m) are no longer eigenvalues of the spin Hamiltonian TL, 
and spin hopping is no longer restricted to ajacent levels. 
Thus transverse field should resolve the bottleneck near 
the top of the barrier, leading to a huge increase of the 
escape prefactor Tq. The aim of the present work is to de- 
scribe this effect by solving the density matrix equation 
(DME) that is a quantum counterpart of the FPE. The 
universal mechanism of spin-lattice relaxation [34, 35] has 
been recently incorporated into the DME [39]. Here it 
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will be used to obtain the results with only one parameter 
describing the spin-phonon interaction, the characteristic 

energy E t = (/wf fi. 3 ) , where p is the mass density of 
the lattice and v t is the speed of transverse sound. The 
model used here is on the same level of simplicity as the 
standard Landau-Lifshitz-Langevin equation or the FPE 
but it is much better justified. We will see that quantum 
effects on the thermal activation rate T do not vanish and 
become even stronger in the large-spin limit S 3$> 1 for 
nearly uniaxial magnetic particles. 

The effective-spin Hamiltonian has the form 



Hs — Ha + Hz, 



(1) 



where Ha is the crystal-field (anisotropy) Hamitonian and 
Hz is the Zeeman Hamiltonian, 
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H z = gm B (H Z S Z + H X S X ) . (2) 



The classical energy barrier AU has particular forms 



AU = DS 2 x 



(l-h x )\ h z = 
{l-h z ) 2 , h x =0, 



(3) 



where h x<z = gm^H Xi Z /(2SD). In general AU(h x ,h z ) can 
be visualized as a Stoner-Wohlfarth astropyramid, com- 
pletely symmetric in h x and h z and basing on the astroid 

In the absence of the transverse field H x , the eigenstatcs 
of the spin are \m) , m = —S, . . . , S, the energy levels being 



e m = —Dm — gm^H z m. 



(4) 



Condition Hu) mm i = e m — e m i = for m ^ ml defines the 
resonance values of the longitudinal field H z : 



gmsH z = kD, 



k = 0,±1,±2, 



(5) 



For these fields all levels in the right well m! = —m — k 
are at resonance with the corresponding levels in the left 
well, m < 0. 

The magnetic particle can be considered as embedded 
in the elastic matrix described by the harmonic-phonon 
Hamiltonian H p h = X)kA ^kA^kA^kA- Approach devel- 
oped in Refs. [34,35] allows to avoid using unknown spin- 
phonon coupling constants and to greatly simplify the for- 
malism. Considering the lattice locally rotated by trans- 
verse phonons without distortion of its crystal field, one 
obtains the spin-phonon interaction 



H s - P h — RHaR 1 — Ha, R — e 



-iS-. 



(6) 



where 5<j) is a small rotation angle given by 5(j) =(1/2) V x 
u(r), u(r) being the lattice displacement due to phonons. 
Expanding Eq. (6) up to first order in 5(f) yields the spin- 
phonon interaction that describes one-phonon processes: 



H 



(i) 



- P h 



H A ,S 



(7) 



It is important that the spin-phonon interaction above 
does not include any poorly known spin-lattice coupling 



coefficients and it is entirely represented by the crystal 
field Ha- To describe the two-phonon (Raman) processes, 
one has to expand i? s - p h up to the second order in 8<j> 
[39,40]. Relaxation rates due to Raman processes are 
generally much smaller than those due to the direct pro- 
cesses since they are the next order in the spin-phonon 
interaction. However, the rates of direct processes can 
be small for special reasons, then Raman processes be- 
come important. Here it happens indeed near the top of 
the barrier in zero transverse field, where the transition 
frequencies between adjacent levels become small. This 
situation has been studied in detail in Ref. [37], however. 
So we will neglect Raman proceses here and concentrate 
on the effect of the transverse field that change transition 
frequencies and drastically increase the escape rate. 

We use the canonical quantization of the lattice dis- 
placement u that yields 



h ^ [ik x e kA ] e ik r / * 

kA 



kA , 



(8) 



Here M is the mass of the unit cell, N is the number of 
cells in the crystal, are unit polarization vectors, A = 
t\,ti,l denotes polarization, and u>k\ = v\k is the phonon 
frequency. Only transverse phonons, ekA-Lk, survive in 
this formula. 

Spin-lattice relaxation including thermal activation can 
be described by the density-matrix equation (DME) [39, 
41]. Early application of the DME to the present model 
in Ref. [42] used the natural basis of states \m) . This pro- 
vided an overall satisfactory description of the thermal 
activation rate, including its strong increase at resonance 
values of H z given by Eq. (5). On the other hand, exact 
energy levels \a) of the spin strongly differ from |m) near 
the top of the barrier even for a small H x . For this rea- 
son, the DME below will be written with respect to the 
energy basis \a) obtained by numerical diagonalization of 
Hs [39]. 

The relaxation terms in the DME can be represented in 
the form that does not explicitly contain Ha, the infor- 
mation about it being absorbed in the spin cigenstates \a) 
and transition frequencies io a p. This can be achieved either 
by changing from the laboratory frame to the local lat- 
tice frame in which Ha remains constant but an effective 
rotation-generated magnetic field arises [26,34,35], or by 
manipulating matrix elements of the spin-phonon interac- 
tion with respect to exact spin states, (a H^ ph f}\ [35]. 

Both methods are mathematically equivalent [35]. As a 
result, the spin part of spin-phonon matrix elements is 
given by the universal expression 



f (i) 

"a/3 



H A ,S 



(a|S|j8) x. g/ i B H. 



p)=iTw a p (a\S\/3) 



(9) 



At tunneling resonances, Eq. (5), one has to use the full 
non-secular form of the DME that couples diagonal ele- 
ments of the density matrix, p aa = n a , to nondiagonal 



p-3 



D. A. Garanin 




Figure 1: Reduced thermal activation rate vs transverse field 
for different temperatures and S = 10, D/ks = 0.65 K. Dra- 
matic increase of F at small h x is due to that of the prefactor To, 
whereas the activation exponent (straight dashed lines) gives 
a comparatively moderate growth. 



elements [39]. In the sequel, tunneling resonances will be 
avoided by choosing the bias field H z in the middle be- 
tween the resonances, to make a better connection with 
classical models. In this case, one can use the system of 
rate equations for the level populations 



dt 



2S+1 

E 

a'=l 



(10) 



where relaxation rates are given by 

/D) 2 [r« (W„'a) + 1) 

+ (aw) n Uai 



(ii) 



Here n u = (e^/C^T) _ -q and 



rW(w) 



0(c), 



(12) 



6(lj) being a Heavyside function and f2 t = (puf/ft) 1 ^ 4 be- 
ing a characteristic frequency. In Eq. (10) transitions 
occur between all the exact spin levels a, although T aa i 
corresponding to pairs of adjacent levels are still domi- 
nating. On the other hand, small transition rates T aa i 
near the top of the barrier are strongly modified even for 
h x <C 1. The coupling of the spin to the environment is 
gauged by a single parameter, £l t in Eq. (12), similarly to 
the paramctrization by the dimcnsionless damping con- 
stant a in the classical LLL equation. However, in the 
present quantum model the rate T^(u>) is frequency de- 
pendent through the distances between energy levels that 
has no analog in the classical scheme. 

Numerical solution of Eq. (10) for the parameters of 
the molecular magnet Mni2 (S — 10, D/k-g, — 0.65 K) 
shifted away from the zero-field resonance, gmsHz = 
0.5D, shows a huge dependence on the transverse field 
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Figure 2: Escape-rate prefactor Fo vs transverse field for dif- 
ferent values of particle's spin S at fixed temperature. 



h x , mainly due to the increase of the prefactor Tq (see 
Fig. 1). The contribution of the Arrhenius exponent 
exp [— AU/ (hsT)] to the growth of T(h x ), shown by 
straight lines exp [2h x DS 2 / (ksT)] following from Eq. (3), 
becomes important only on the right side of the plot where 
the growth of To(h x ) saturates. The effect of the trans- 
verse field here is much greater than in the classical model, 
the LD curve in Fig. 3 of Ref. [38]. Note that in the present 
model we are in the uniaxial - low damping limit since the 
damping calculated here from the first principles for real- 
istic f2 t is much smaller than all other frequency scales. 

For effective spins of magnetic particles that are much 
greater than S — 10, the effect of the transverse field is 
huge. Since T in zero transverse field becomes anomalously 
small for large spins, one cannot normalize the results by 
it. It is better to plot the prefactor To alone defined as 
To = r exp [AU I (hsT)] , where T follows from the solution 
of Eq. (10) and AU is found numerically for the classical 
model. The characteristic rate 



r = sr s .-i 



J S,S-1 



12ttO| 



(13) 



can be used to normalize the results for Tq in a wide 
range of h x . Here Ts,s—i an d ws,S-i are zero-temperature 
relaxation rate and transition frequency for the lowest- 
lying pair of levels in the well, defined above. T = 
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STs,s—i is an overall measure of the relaxation rate in- 
side a well. Indeed, in the natural basis the spin-phonon 
transition rate between two adacent levels is proportional 



t0 4,m±i I 42 ], where I 



m,m±l 



'm,m±i(2m ± 1) and 



,JS(S + 1) -m{m±l). In r S:S _i the factor 
i m m±1 yields S while for a typical value of m in the in- 



^m.m±l 

terval — S < m < S it yields S 2 . This is the origin of an 
additional S in Eq. (13). 

In the comparizon between different values of S shown 
in Fig. 2 the product SD is kept constant, as it should 
be for the effective anisotropy of magnetic particles. In 
numerical calculations D/k-Q = 6.5/ S K is used and the 
temperature k B T = SD = const. Fig. 2 (a) shows that 
curves r /r scale for large S in a wide range of h x , that 
means T ct S 2 . Calculations use custom- precision matrix 
algebra within Wolfram Mathematica and become slow for 
spins as large as S — 80. One can see that in the large- 
spin limit To becomes small if h x — > and h x — * 1- In 
particular, for h x — *■ the apparent behavior is Tq cc h 2 . 

The behavior of Tq at small transverse fields is eluci- 
dated in Fig. 2 (b). Here one has to use a slightly dif- 
ferent normalization of Tq to make curves collapse in a 
wide range of h x , yielding r oc S 3 / 2 h 2 . In the uniaxial 
limit h x — > the curves for different S diverge. Here the 
escape prefactor is given by the transition rate between 
the ajacent levels near the top of the barrier r mjTO ±i with 
m ~ 1. Using Eq. (A9) of Ref. [35] for T„ um ±i [multiplied 
by n Umm ± 1 = k&T / (fou} m m ±i) to account for a nonzero 
temperature], one obtains Tq cc S~ 2 . This is the top-of- 
the-barrier bottleneck mentioned in the introduction. In 
the representation of Fig. 2 (b) one has Tq/S 3 / 2 cc S~ 7 ^ 2 . 
One can see that doubling S results in the drop by a factor 
2" 7 1 2 ~ 11 in the asymptotic h x — > values in Fig. 2 (b). 

The anomalously small rate in the uniaxial limit above 
is in part due to the factor 2m±l discussed below Eq. (13). 
In the zero-bias case the top of the barrier corresponds to 
m ~ 1 that results in additional smallncss. In the case of 
a strong enough bias one has 2m ± 1 ~ S near the top 
of the barrier, so that the anomalously small escape rate 
solely results from small Lu m .m±i- The results of numerical 
calculations for the bias h x w 0.5, adjusted to the middle 
between two tunneling resonances are shown in Fig. 3. 
The curves for r in a broad range of h x in Fig. 3 (a) look 
complicated for moderate S but still collapse for large S. 
The decrease of Tq at h x — ► is indeed weaker than in the 
unbiased case above. The results at small h x in Fig. 3 (b) 
show a dependence Tq cc S 3 / 2 h®- S5 : where the exponent 
0.85 cannot be easily explained. For h x = Eq. (A9) 
of Ref. [35] yields Tq cc S° in the biased case, also much 
smaller than Tq cc S 2 for h x ~ 1. 

It should be noted that the secular approximation lead- 
ing to Eq. (10) relies on the smallness of the relaxation 
terms in the DME in comparison to the dissipationless 
terms for the nondiagonal elements of the density ma- 
trix [39] . Then slow diagonal elements p aa — n a dynam- 
ically decouple from the fast nondiagonal terms p a p. The 
classical match of the secular approximation is the low- 
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Figure 3: Escape-rate prefactor To vs transverse field for dif- 
ferent values of particle's spin S at fixed temperature at strong 
bias, h x « 0.5 (between tunneling resonances). 



damping (LD) approximation introduced by Kramers for 
a particle in a potential well [43]. In the LD limit the en- 
ergy of the particle or spin is nearly conserved, so that the 
fast motion over constant-energy trajectories averages out 
and what is left is the slow energy diffusion (see, e.g., Eqs. 
(15) and (16) of Ref. [38]). Similarly, Eq. (10) describe a 
slow hopping over the quantum energy levels of the spin. 

One can ask whether the richness of damping regimes 
that exist in the classical-spin model [38] can be realized 
for a realistic large quantum spin of a magnetic particle. 
For instance, the intermediate-to-large damping (IHD) 
case requires that the gyroscopic and relaxation terms in 
the FPE be comparable. This means that the dissipation- 
less and relaxation terms in the DME be comparable as 
well. Of course, for S ^> 1 nondiagonal elements p a p close 
to diagonal become slow as ui a f3- However, the relaxation 
rate T a p between the states a and /3 scales as T a p cc ui 2 a Q 
for \u) a p\ -C T and decreases faster than uj a j3 in the quasi- 
classical limit S>1. 

What can change the situation is Raman processes that 
become independent of u) a p for small u a p. Incorporating 
Raman processes requires generalization of the results of 
Ref. [37] for a nonzero transverse field that is a nontrivial 
task. As Raman processes are much weaker than direct 
processes, crossover to a Raman-dominated behavior re- 
quires very large S. Although a general nonsecular DME 
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can be solved as described in Ref. [39], calculations arc 
much slower than those of Eq. (10) and become prohibitive 
for the required very large S. For this reason, Raman pro- 
cesses cannot be adequately treated within this Letter and 
should be considered elsewhere. 

In all cases, even with account of Raman processes, 
there should be a bottleneck for spin diffusion near the 
top of the barrier in the case of nearly uniaxial magnetic 
particles. Transverse magnetic field gradually resolves the 
bottleneck and leads to a huge increase of the escape-rate 
prefactor Tq that is more important than the barrier low- 
ering. This is the main finding of this work. Thermal 
activation rates of nearly-uniaxial magnetic particles are 
very sensitive to any deviations from the axial magnetic 
symmetry, e.g., due to surface anisotropy [15-17]. Robust 
results require a strong enough transverse field. 

Whether for large spins the classical stochastic ap- 
proach could be modified to embrace the spacing between 
quantum-mechanical levels that has been shown to be im- 
portant, remains an open question. 

This work has been supported by the Cottrell College 
Science Award of the Research Corporation. 
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